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I. INTRODUCTION 

In this chapter we will discuss an application of pre- 
cision molecular spectroscopy to the studies of the pos- 
sible spatial and temporal variations of the fundamen- 
tal constants. As we will see below, molecular spec- 
tra are mostly sensitive to two such dimensionless con- 



stants, namely the fine-structure constant a 



h c 



and 



the electron-to-proton mass ratio /i = m e /m p (note that 
some papers define [i as an inverse value, i.e. proton- 
to-electron mass ratio). At present NIST gives following 
values of these constants pQ: a' 1 = 137.035999679(94) 
and pT l = 1836.15267247(80). 

The fine-structure constant a determines the strength 
of electromagnetic (and more generally electroweak) in- 
teractions. In principle, there is similar coupling constant 
a s for quantum chromodynamics (QCD). However, be- 
cause of the highly nonlinear character of the strong in- 
teractions, this constant is not well defined. Therefore, 
the strength of the strong interactions is usually charac- 
terized by the parameter Aqcd , which has the dimension 
of mass and is defined as the position of the Landau pole 
in the logarithm for the running strong coupling constant, 
a s (r) = const/ In (r Aqcd /he), where r has dimension of 
length. 

In the Standard Model (SM) there is another fun- 
damental dimensional parameter — the Higgs vacuum 
expectation value (VEV), which determines electroweak 
unification scale. Electron mass m e and quark masses 
m q are proportional to the Higgs VEV. Consequently, 
the dimensionless parameters X e = m e / 'Aqcd and X q — 
to 9 /Aqcd link electroweak unification scale with strong 
scale. For the light quarks u and d, X q <C 1. Because 
of that the proton mass m p is proportional to Aqcd and 
X e oc /i. Below we will use /i instead of X e because it is 
more directly linked to experimentally measured atomic 
and molecular observables. 

Below we will show that huge enhancement of the 
relative variation happens in transitions between close 
atomic, molecular and nuclear energy levels. Recently 
several new cases were found, where the levels are very 
close and narrow. Large enhancement of the variation 
effects is also possible in cold collisions of atoms and 
molecules near Feshbach resonances. 

We will start with general review of the present situ- 
ation in the search of the variation of a and After 
that we will discuss in more detail the results, which 
follow from the astrophysical observations of the opti- 



cal and microwave spectra of molecules. Finally, we will 
describe possible laboratory experiments with molecules. 
This field is very new and there are no competitive labo- 
ratory results on time- variation with molecules yet (see, 
however, Sec. VII I, but there are very promising propos- 



als and several groups already started experiments. 

The analysis of the data from Big Bang nucleosynthesis 
PJ, quasar absorption spectra, and Oklo natural nuclear 
reactor give us the space-time variation of constants on 
the Universe lifetime scale, i.e. on times from few bil- 
lion to more than ten billion years. Comparison of the 
frequencies of different atomic and molecular transitions 
in laboratory experiments gives us the present variation 
on the timescale from few months to few years. There is 
no model independent connection between variations on 
such different timescales. However, in order to compare 
the importance of different results, we will often assume 
linear time dependence of the constants. This way we 
can interpret all results in terms of time derivatives of 
the fundamental constants. Within this assumption, we 
can use quasar absorption spectra to obtain the best cur- 
rent limit on the variation of the mass ratio /i and X e [3] : 



Mli = X e /X e = (1 ± 3) x 10- 16 yr- 1 



(1) 



We can also combine this result with the atomic clock 
results [J to get the best limit on variation of a [5J El H] : 



a/a = (-0.8 ± 0.8) x 10~ 16 yr" 1 



(2) 



Note that both limits ([lj and (2]) depend on the assump- 
tion of the linear time dependence of fundamental con- 
stants. 

The Oklo natural reactor gives the best limit on the 
variation of X s — m s /AQCD where m s is the strange 
quark mass [SI El HD]: 



\X S /X S 



< 10~ 18 yr" 1 . 



(3) 



Note that the Oklo data can not give us any limit on 
the variation of a since the effect of a there is much 
smaller than the effect of X s and within the accuracy of 
the present theory should be neglected [10] . 

In addition to the time-variation, one can also con- 
sider spatial- variation of constants. Massive bodies (stars 
or galaxies) can also affect physical constants. In other 
words the fundamental constants may depend on the 



2 



gravitational potential, e.g. 

da/a = k a S(GM/rc 2 ) , (4) 

where G is the gravitational constant and r is the distance 
from the mass M. The strongest limit on such variation: 

fc Q + 0.17fc M = (-3.5 ±6) x 1CT 7 , (5) 

is obtained in Ref. [6] from the measurements of the de- 
pendence of atomic frequencies on the distance from the 
Sun due to the ellipticity of the Earth's orbit [HE] (pa- 
rameters is defined by analogy with Eq. Q). Be- 
low we will also discuss some other results, including 
those, which indicate nonzero variation of fundamental 
constants. 



II. THEORETICAL MOTIVATION 

How changing physical constants and violation of local 
position invariance may occur? Light scalar fields very 
naturally appear in modern cosmological models, affect- 
ing parameters of SM including a and /i (for the whole list 
of SM parameters see [H]). Cosmological variations of 
these scalar fields should occur because of drastic changes 
of the composition of the Universe during its evolution. 

Theories unifying gravity and other interactions sug- 
gest the possibility of spatial and temporal variation of 
physical "constants" in the Universe [13]. Moreover, 
there exists a mechanism for making all coupling con- 
stants and masses of elementary particles both space and 
time dependent, and influenced by local environment (see 
review [14). Variation of coupling constants can be non- 
monotonic, such as damped oscillations, for instance. 

These variations are usually associated with the effect 
of massless (or very light) scalar fields. One candidate 
is the dilaton: a scalar which appears in string theories 
together with graviton, in a massless multiplct of closed 
string excitations. Other scalars naturally appear in cos- 
mological models, in which our Universe is a "brane" 
floating in a space of larger dimensions. The scalars are 
simply brane coordinates in extra dimensions. However, 
the only relevant scalar field recently discovered, the cos- 
mological dark energy, so far does not show visible vari- 
ations. Observational limits on the variations of physical 
constant given in Sec. [I] are quite stringent, allowing only 
scalar couplings, which are tiny in comparison with grav- 
ity. 

A possible explanation was suggested by Damour et 
al. [T5J [TH] who pointed out that cosmological evolu- 
tion of scalars naturally leads to their self-decoupling. 
Damour and Polyakov have further suggested that varia- 
tions should happen when the scalars get excited by some 
physical change in the Universe, such as phase transi- 
tions, or other drastic changes in the equation of state of 
the Universe. They considered several of them, but since 
the time of their paper a new fascinating transition has 
been discovered: from matter dominated (decelerating) 



era to dark-energy dominated (accelerating) era. It is 
a relatively recent event, corresponding to cosmological 
rcdshift z w 0.5, or the look-back time of approximately 
5 billion years. 

The time dependence of the perturbation related to 
this transition can be calculated, and it turned out 
[i"7] li"8] that the self-decoupling process is effective 
enough to explain why after this transition the varia- 
tion of constants is as small as observed in laboratory 
experiments at the present time, while being at the same 
time consistent with possible observations of the vari- 
ations of the electromagnetic fine-structure constant at 

z > l [HJ[2D1[2T]. 



III. DEPENDENCE OF ATOMIC AND 
MOLECULAR SPECTRA ON a AND /x 

Atomic and molecular spectra are most naturally de- 
scribed in atomic units {Ti — m e = e = 1), where energy 
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is measured in Hartrees (1 Hartree = e ^" e = 2 Ry = 
219474.6313705(15) cm -1 ). In these units nonrelativistic 
Schrodinger equation for an atom with infinitely heavy 
pointlike nucleus does not include any dimensional pa- 
rameters. The dependence of the spectrum on a appears 
only through relativistic corrections, which describe fine- 
structure, Lamb shift, etc. The dependence of atomic 
energies on /j, is known as isotope effect and is caused by 
finite nuclear mass and volume. There are even smaller 
corrections to atomic energies, which depend on both a 
and /i and are known as hyperfine structure. 

One can argue that atomic energy unit itself depends 
on a as it can be expressed as a 2 m e c 2 , where m e c 2 is 
the rest energy of the free electron. However, experimen- 
tal search for possible variation of fundamental constants 
consists in observing time-variations of the ratios of dif- 
ferent transition frequencies to each other. In such ratios 
the dependence of the units on fundamental constants 
cancels out. Below we will use atomic units unless oth- 
erwise is explicitly stated. 

Relativistic corrections to the binding energies of 
atomic valence electrons are of the order of a 2 Z 2 , where 
Z is atomic number and become quite large for heavy 
elements. For our purposes, it is convenient to present 
the dependence of atomic transition frequencies on a 2 in 
the form 

u) = ujo + qx, (6) 

where x — ( — ) 2 — 1 « and wn is a transition fre- 
quency for a = ao- Rough estimates of q- factors can be 
obtained from simple one-particle models, but in order to 
obtain accurate values one has to account for electronic 
correlations and perform large-scale numerical calcula- 
tions. Recently such calculations were made for many 
atoms and ions [H] OH M HH1 123 OH OH OS] ■ 

Isotope effects in atoms are of the order of fi ~ 
10~ 3 and magnetic hyperfine structure roughly scales as 
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a 2 [iZg nuc ~ 10~ 7 Zg nuc , where g nuc is nuclear g-factor. 
One has to keep in mind that g nuc also depends on \i and 
quark parameters X q . This dependence has to be consid- 
ered, when we compare, for example, the frequency of the 
hyperfme transition in 133 Cs (Cs frequency standard) [5J, 
or the hydrogen 21 cm hypcrfinc line [30 , 31j to various 
optical transitions [5J. 

At present there are many very accurate experiments 
where different optical and microwave atomic clocks are 
compared to each other [H E2 E3 EH 1311 E3 EH EH ES] . 
These experiments place strong limits on the time- 
variation of different combinations of a, fi, and g n uc- As 
we mentioned above, the limit on a-variation Q follows 
from the experiment [4] and the limit (JT|) in the assump- 
tion of linear time-dependence of all constants. A de- 
tailed discussion of atomic experiments can be found in 
recent reviews [40l 141 j . 

On a cosmological timescale a comparison of the hy- 
pcrfinc transition in atomic hydrogen with optical tran- 
sitions in ions, was done in Refs. [30] EI]- This method 
allows one to study time-variation of the parameter 
F = a 2 g p ^jL, where g p is proton g-factor. Analysis of 
the absorbtion spectra of nine quasars with redshifts 
0.23 < z < 2.35 gave 

6F/F= (6.3 ±9.9) x 10~ 6 , (7) 
F/F= (-6± 12) x 10~ 16 yr" 1 , (8) 

which is consistent with zero variation of /z and a. 

Molecular spectroscopy opens additional possibilities 
to study variation of fundamental constants. It is known 
that /i defines the scales of electronic, vibrational, and ro- 
tational intervals in molecular spectra, E e \ : E v ^ : E TOt ~ 
1 : /i 1 / 2 : /i. In addition to that, molecules also have fine 
and hyperfme structure, A-doubling, hindered rotations, 
etc. All these structures have different dependencies on 
fundamental constants. Obviously, comparison of these 
structures to each other allows the study of different com- 
binations of fundamental constants. 

The sensitivity to temporal variation of fundamental 
constants may be strongly enhanced in transitions be- 
tween narrow close levels of different nature. Huge en- 
hancement of the relative variation Slo/lo can be obtained 
in transition between almost degenerate levels in atoms 
HH M H3 HI S3], molecules [3J Si S3 Si SO], and 
nuclei [501 ED- 

An interesting case of enhancement of the effect of vari- 
ation of fundamental constants can be found in collisions 
of ultracold atoms and molecules near Feshbach reso- 
nances [55]. The scattering length A near the resonance 
is extremely sensitive to the /i- variation: 



SA 5fi 

= & — 

A fl 



(9) 



where the enhancement factor K can be very large. For 
example, for Cs-Cs collisions K ~ 400 [55]. Enhancement 
can be further increased by adjusting the position of the 



resonance using external fields. Near a narrow magnetic 
or optical Feshbach resonance the enhancement factor K 
may be increased by many orders of magnitude. 

Calculation of the factor K in Ref. [55] is based on 
the analytical formula for the scattering length derived 
in Ref. |53j . This formula is valid for an arbitrary in- 
teratomic potential with power long-range tale (l/r 6 for 
neutral atoms), i.e. this result includes all unharmonic 
corrections. 

To the best of our knowledge, it is the only suggested 
experiment on time- variation, where the observable is not 
frequency. Because of that, we have to find another pa- 
rameter L of the dimension of length to compare A with. 
In Ref. [55] the scattering length was defined in atomic 
units (as). It is important, however, that because of the 
large enhancement in Eq. ([9]), the possible dependence of 
L on fi becomes irrelevant. For example, if we measure 
A in conventional units, meters, which are linked to Cs 
standard, then SL/L = —Sfi/fi, and 



(A/L) [i 



(10) 



As long as K ^> 1 the dependence of the used units on 
fundamental constants can be neglected. Below we will 
discuss several other experiments with huge enhancement 
factors, where this argument can be also applied. 



IV. ASTROPHYSICAL OBSERVATIONS OF 
THE SPECTRUM OF H 2 



H2 is the most common molecule in the Universe and 
its UV spectra have been used for the studies of the pos- 
sible /i- variation for a long time. For a given electronic 
transition, the frequency of each rovibrational line has 
different dependence on /j, [5H[55]. Therefore, compar- 
ison of rovibrational frequencies from astrophysics with 
laboratory observations can give information on fi. 

In the adiabatic approximation, the rovibrational lev- 
els of the electronic state A with vibrational and rota- 
tional quantum numbers v and J are given by the Dun- 
ham expression [55]: 



E(v,J)= J2 + 



A 



(11) 



fe,/>0 



where each term depends on /1 in a following way: 

Y ktl cx n l+k / 2 . (12) 

Because of the smallncss of the parameter /1, coefficients 
Y^ i rapidly decrease with both k and Z, and for small v 
and J, we have the usual vibrational (k = 1) and rota- 
tional (I = 1) terms. The zero term of this expansion 
(k = I = 0) corresponds to the electronic energy. 

One can define the sensitivity coefficient K t for each 
rovibrational transition i of a given electronic band e — g 
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as we have seen, the sensitivity coefficients Ki in Eq. ( 13 I 



dvj \ / / dp 
n J I \ /x 
/x / (i-Ee di? s 



E K - E. 



(13) 



where both energies are given by expansion (111. The 



sign of Ki depends on the rovibrational energies of the 
excited (e) and ground (g) states (in the absorbtion spec- 
tra of the quasars only transitions from the ground elec- 
tronic state are seen). The electronic energy, presented 
by the term Y 00l dominates the expansion and the co- 
efficients Ki are rather small. Typically they are on the 
order 1CP 2 , but can reach 0.05 for large values of the 
quantum numbers v and J. 



The coefficients of expansion (11) can be found by fit- 
ting experimental spectra. After that the sensitivity co- 
efficients Ki are found from Eqs. (12 1 and (13). Some 



rovibrational levels of different electronic excited states 
appear to be very close. For such levels additional non- 
adiabatic corrections can be included within the two-level 
approximation [57] , 

If there is /i-variation Afi, this would lead to a differ- 
ence in observed redshifts z, for different lines: 



G 



^q.abs A/2 ^ 



1 ~h ^q,abs 



(14) 



By plotting the reduced redshifts Q against the sen- 
sitivity coefficients Ki, one can estimate A////J. The 
most recent study [20 of the possible /x-variation us- 
ing astrophysical data on H2 was based on the observa- 
tion of the two quasar absorbtion systems with redshifts 
2 q ,abs = 3.02 and 2.59. An analysis of the data on 76 
lines from two UV bands of H2 gave the following result: 



= (-20 ±6) x 10" 



(15) 



This result indicates, at a 3.5(7 confidence level, that /i 
has increased during the past 12 billion years. Assuming 
linear time-dependence we can rewrite Eq. (15 1 as 



/' 



= (17 ±5) x 10 



-16 



vi- 



ae) 



This has to be compared with the ammonia result (JT|), 
which corresponds to a timescale about 6.5 billion years 
and is discussed in more detail in Sec. IVI1 



are rather small. This is because of the relative smallness 
of rovibrational energy compared to the total transition 
energy. Thus, it may be useful to study microwave spec- 
tra of molecules, where the relative frequency variations 
due to varying constants are larger. 



A. Rotational spectra 

In 1996 Varshalovich and Potekhin [58 compared red- 
shifts for microwave rotational transitions (J = 3 — > 
J = 2) and (J = 2 -> J = 1) in the CO molecule 
with redshifts of optical lines of light atomic ions from 
the same astrophysical objects at redshifts z — 2.286 
and z — 1.944. As long as atomic frequencies are inde- 
pendent on /i and rotational transition frequencies arc 
proportional to /i, this comparison allowed to set the fol- 
lowing limits on variation of //: 



A 4 
H 



= (-0.6 ±3.7) x 10~ 4 at z = 2.286, (17a) 
= (-0.7 ± 1.0) x 10" 4 at z = 1.944. (17b) 



In the same paper |58j . the authors compared the 
(J = — ► J = 1) CO absorbtion line with the 21 cm 
hydrogen line for an object with z = 0.2467. They did 
not find a significant difference in respective redshifts and 
interpreted this result as yet another limit on variation 
of (i. However, as we mentioned above, the frequency 
of the hydrogen hypcrfine line is proportional to a 2 [ig pi 
and this result actually places limit on the variation of 
the parameter F — a 2 g p [59 . Recently a similar analysis 
was performed by Murphy et al. [60 using more accurate 
data for the same object at Z = 0.247 and for a more dis- 
tant object at z = 0.6847, and the following limits were 
obtained: 



5F 



(-2.0 ±4.4) x 10" 
(-1.6 ±5.4) x 10" 



at z = 0.2467 , (18a) 
at z = 0.6847 . (18b) 



The object at z = 0.6847 is associated with the gravita- 
tional lens toward quasar B0218+357 and corresponds to 
the backward time ~ 6.5 Gyr. This object was also used 
by other authors, as will be discussed in Sec. |VB| and 
Sec. IVTl 



V. ASTROPHYSICAL OBSERVATIONS OF 
MICROWAVE MOLECULAR SPECTRA 

In the previous section we discussed astrophysical ob- 
servations of UV spectra of H 2 . The corresponding ab- 
sorbtion bands are very strong and can be observed even 
for objects with very high redshifts. On the other hand, 



B. The 18 cm transitions in OH 

Let us consider transitions between hyperfine substates 
of the 2 n 3 / 2 ground-state A-doublet in the OH molecule 
(61] [62l [63] . The A-doubling for 2 n 3 / 2 states appears in 
the third order in the Coriolis interaction and the corre- 
sponding energy interval is inversely proportional to the 
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spin-orbit splitting between the 2 n 3 / 2 and 2 n 1 / 2 states, 
i.e., it scales as /j, 3 a~ 2 , while the hyperfine structure in- 
tervals scale as a 2 [ig nuc . Therefore, the ratio of the hy- 
perfine interval to the A-doubling interval depends on the 
combination F = a 4 /i~ 2 <? n uc- Higher-order corrections 



modify this parameter to the form F 



,,3.14 



-1.57 



f-Ju 



The hyperfine-structure splitting for the OH molecule 
is approximately 50 MHz and is much smaller than A- 
doubling interval, which is about 1700 MHz. Because of 
that, it is actually easier to compare the A-doubling tran- 
sitions in OH to the 21 cm hydrogen line, or to rotational 
lines of the HCO+ molecule [§3 El M IM] • 

The most stringent limit on the variation of F was 
obtained in Ref. [M] from observations of the z — 0.6847 
gravitational lens: 



AF/F = (0.44 ± 0.36 stat ± 1.0 syst ) x 10" 



(19) 



where systematic error mostly accounts for the possible 
Doppler noise, i.e. for the possible difference in the ve- 
locity distributions of different molecules in a molecular 
cloud. 

The laboratory frequencies of the OH A-doublet were 
recently rcmcasurcd with higher precision using cold 
molecules produced by a Stark decelerator [55] ■ That 
may become important for future astrophysical measure- 
ments with higher accuracy. 



VI. LIMIT ON TIME- VARIATION OF fi FROM 
INVERSION SPECTRUM OF AMMONIA 

Several years ago, van Veldhoven et at. suggested to 
use decelerated molecular beam of ND3 to search for the 
variation of fi in laboratory experiments [46 . The am- 
monia molecule has a pyramidal shape and the inversion 
frequency depends on the exponentially small probability 
of tunneling of the three hydrogen (or deuterium) atoms 
through the potential barrier [55] . Because of that, it is 
very sensitive to any changes of the parameters of the 
system, particularly to the reduced mass for this vibra- 
tional mode. The authors of Ref. [46 found that for 
ND3 molecule, Suj/lu = — 5.6<5/i//i. Therefore, the inver- 
sion frequency of ND3 is nearly an order of magnitude 
more sensitive to /x-variation than typical molecular vi- 
brational frequencies (note that Ref. [46 contains a mis- 
print in the sign of the effect). 

However, even such enhanced sensitivity is insufficient 
to make the laboratory experiment on the time- variation 
of [i using conventional molecular beams competitive. 
Stark-deceleration was used in Ref. [IB] to slow down the 
beam to 52 m/s. Still, a much slower beam, or a fountain 
is necessary to increase the sensitivity by several orders of 
magnitude before a competitive experiment can be per- 
formed. The work in this direction is in progress |67j . 

On the other hand, an only slightly smaller enhance- 
ment also exists for the inversion spectrum of NH3, which 



is often seen in astrophysics, even for high z objects. This 
fact was used in [3] to place the limit ([lj, which we will 
now discuss in some detail. The inversion vibrational 
mode of ammonia is described by a double- well potential 
with the first two vibrational levels lying below the bar- 
rier. Because of the tunneling, these two levels are split in 
inversion doublets. The lower doublet corresponds to the 
wavelength A w 1.25 cm and is used in ammonia masers. 
Molecular rotation leads to the centrifugal distortion of 
the potential curve. Because of that, the inversion split- 
ting depends on the rotational angular momentum J and 
its projection K on the molecular symmetry axis: 



Ci [j{J+\)-K 2 ] +c 2 K 2 



(20) 



where we omitted terms with higher powers of J and K. 
Numerically, d? nv w 23.787 GHz, a w 151.3 MHz, and 
c 2 « 59.7 MHz. 



In addition to the rotational structure ( 20 ) the inver 



sion spectrum includes much smaller hyperfine structure. 
For the main nitrogen isotope 14 N, the hyperfine struc- 
ture is dominated by the electric quadrupole interaction 
(~ 1 MHz) [68]. Because of the dipolc selection rule 
AK = the levels with J = K are metastable and 
in laboratory beam experiments the width of the cor- 
responding inversion lines is usually determined by colli- 
sional broadening. In astrophysics the lines with J = K 
are also narrower and stronger than others, but the hy- 
perfine structure for spectra with high rcdshifts is still 
unresolved. 

For our purposes it is important to know how the pa- 



rameters in Eq. (20) depend on fundamental constants. 



The molecular electrostatic potential in atomic units does 
not depend on fundamental constants (here we neglect 
small rclativistic corrections which give a weak a depen- 
dence). Therefore, the inversion frequency u>f nv and con- 
stants C12 are functions of [i only. Note that the coeffi- 
cients Ci depend on /i through the reduced mass of the 
inversion mode and because they are inversely propor- 
tional to the molecular moments of inertia. This implies 
a different scaling of tu? nv and Cj with /i. 



The inversion spectrum ( 20 1 can be approximately de 



scribed by the following Hamiltonian: 



FliyiV 



2Afi u x 



U(x) 



(21) 



[J{J- 



1) 



K 2 



I 2 (x) 



K 



where x is the distance from N to the H-plane, / l5 I 2 
are moments of inertia perpendicular and parallel to the 
molecular axis, correspondingly, and M\ is the reduced 
mass for the inversion mode. If we assume that the length 
d of the N — H bond does not change during inversion, 
then Mi — 2.54 m p and 



l2{x) » 3m p d 2 [l 



\m v d 2 [1 + 0.2(x/d) 2 ] 
(x/d) 2 } . 



(22) 
(23) 
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The dependence of l\ t % on x generates correction to the 
potential energy of the form C(J, K) x 2 fx. This changes 
the vibrational frequency and the effective height of the 
potential barrier, therefore changing the inversion fre- 
quency tt>inv given by Eq. ( 20 1 . 



get: 



SwL 



6 ft 
2^ 



l + S- 



S W v ib 



2 U B 



E 



4.4 



(29) 



which is close to the numerical result ( 25 I 



Following Rcf. (69] we can write the potential U(x) in 
Eq. (21) in the following form: 



U(x) 



2 kx 



&exp (- 



(24) 



We see that the inversion frequency of NH3 is an order 
of magnitude more sensitive to the change of fj, than typ- 
ical vibrational frequencies. The reason for this is clear 
from Eq. ( 29 ) : it is the large value of the action S for the 



Fitting vibrational frequencies for NH3 and ND3 gives 
k m 0.7598 a.u., b « 0.05684 a.u., and c » 1.3696 a.u. 
Numerical integration of the Schrodinger equation with 
the potential (24 1 for different values of fi gives the fol- 
lowing result: 



tunneling process. 



Using Eqs. (21) - ( 23 > one can also find the dependence 
on [i of the constants Ci,2 in Eq. ( p0| [3J: 



C1.2 



5.1 



Sfi. 



(30) 



4.46 



5fi 



(25) 



It is instructive to reproduce this result from an analyti- 
cal calculation. In the WKB approximation the inversion 
frequency is estimated as [70] : 



It is clear that the above consideration is directly ap- 
plicable to ND3, where the inversion frequency is 15 
times smaller and Eq. (26b I gives S « 8.4. According to 
Eq. ( 29 1 , this leads to a somewhat higher sensitivity of 



the inversion frequency to \i in agreement with Ref. 



^vib 

7T 
^vib 

7T 



exp (-S) 



(26a) 



ND, 



exp 



a/2Mi(E/(z) - E) dx 



Sc 2 

C2 



5.7 



6.2 



(31) 



(26b) 



We see from Eqs. (251 and (30 1 that the inversion fre- 



quency cj; nv and the rotational intervals 



where w V ib is the vibrational frequency of the inversion 
mode, S is the action in units of K, x = ±<z are classical 
turning points for the energy E. For the lowest vibra- 
tional state E = Umin - 
values w v ib = 950 cm" 
S&5.9. 



Wi nv ( J2, K2) have different dependencies on /i. In princi- 
ple, this allows one to study time-variation of fi by com- 
paring different intervals in the inversion spectrum of am- 
monia. For example, if we compare the rotational interval 



^Wvib- Using the experimental 

and uj inv = 0.8 cm -1 , we get to the inversion frequency, then Eqs. |25J and (|30j) give 



Expression (26b) allows one to calculate the depen- 
dence of u>f nv on the mass ratio //. Let us present S in the 
following form: S — A^ 1 ! 2 J" \fU{x) — Edx, where A 
is a numerical constant and the square root depends on 
fi via E: 



^inv 

[uii nv (Ji, Ki) 

^inv (J 2 ,if 2 )]/< v 



0.6 



(32) 



1 

2/t 
1 

2fi 



ds; 
as 



d^dE 
dE~dii 



(27a) 
(27b) 



The relative effects are substantially larger if we compare 
the inversion transitions with the transitions between the 
quadrupole and magnetic hyperfine components. How- 
ever, in practice, this method will not work because of 
the smallness of the hyperfine structure compared to typ- 
ical linewidths in astrophysics. 

Again, as in the case of A-doubling in OH molecule, it 
is more promising to compare the inversion spectrum of 
NH 3 with rotational spectra of other molecules, where 



It is easy to see that dS/dfi = —S/2[i. The value of 
the third term in Eq. ( 27b ) depends on the form of the 
potential barrier: 



5u> T , 



Wrot 



A 4 



(33) 



as 

dE 



s 



In astrophysics any frequency shift is related to a corre- 
sponding apparent redshift: 



4(7„ 



E 



(28) 



where for a square barrier q = 1, and for a triangular 
barrier q — 3. For a more realistic barrier shape, g«2. 
Using the parametrization ( 24 1 to determine U max we 



5u> 



l + z 



(34) 



According to Eqs. (25 1 and (33 1, for a given astrophysical 
object with z = zq variation of fi leads to a change of the 
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apparent redshifts of all rotational lines Sz YO t — — (1 + 
zq) Sfi/ fj, and the corresponding shifts of all inversion lines 
of ammonia are: 8z- lm — —4.46 (1 + z Q ) Sfi/fi. Therefore, 
comparing the apparent redshift z- nlv for NH3 with the 
apparent redshifts z rot for rotational lines we can find 
5fi/fi: 



/' 



0.289 



Zrot 



1 + ^0 



(35) 



High-precision data on the redshifts of NH3 inversion 
lines exist for the already mentioned object B0218+357 
at z ss 0.6847 [71]. Comparing them with the redshifts of 
rotational lines of CO, HCO + , and HCN molecules from 
Rcf. \72\ one can get the following conservative limit from 
Eq. (3l: 



= (-0.6 ± 1.9) x 10" 



(36) 



Taking into account that the redshift z w 0.68 for the 
object B02 18+357 corresponds to the look-back time of 
about 6.5 Gyr, this limit translates into the most strin- 
gent present limit (JTJ for the variation rate \ij \i. 



([IJ, but there are good chances that it will be soon sig- 
nificantly improved. 



VIII. CLOSE NARROW LEVELS IN DIATOMIC 
MOLECULES 



In this section we focus on very close narrow levels of 
different nature in diatomic molecules. Such levels may 
occur due to cancelation between either hyperfine and 
rotational structures [45], or between the fine and vibra- 
tional structures of the electronic ground state [49] . The 
intervals between the levels correspond to microwave fre- 
quency range convenient for experiments and the level 
widths are very small, typically ~ 10 -2 Hz. The en- 
hancement of the relative variation K can exceed 10 5 . 



Molecules with cancelation between hyperfine 
structure and rotational intervals 



VII. EXPERIMENT WITH SF 6 

Now we switch to laboratory molecular experiments on 
time-variation. We start with the recent experiment on 
two-photon vibrational transition (v = 0, J = 4) — ► (v = 
2, J = 3) in SFg 07] . This is a Ramsey-type experiment 
with a supersonic beam of SFg molecules. The beam 
velocity u — 400 m/s and the length of the interaction 
region D = 1 m corresponds to the linewidth of uj ID — 
200 Hz. 

A CO2 laser was used to drive the two-photon tran- 
sition and its frequency was controlled by a Cs stan- 
dard [?5]. This means, that the vibrational frequency 
tj V ib in SFg was compared with the hyperfine transition 
frequency Whfs in Cs. Therefore, the experiment was 
sensitive to the combination of fundamental constants 



F 



-l/2 a 2.83^ 



Measurements continued for 18 



months, and the following result was obtained: 



F/F 



(1.4 ±3.2) x lO^yr" 1 



(37) 



This limit is weaker, than the most stringent limits 
obtained with atomic clocks. On the other hand, it con- 
strains a different combination of fundamental parame- 
ters. Most importantly, in atomic experiments the pa- 
rameters g n and /i always go as a product while 
here we have combination g n [i~ x l 2 . That allows to com- 



bine atomic results [4] [35] [37] with limit (37 1 to obtain 
the best laboratory limit on /x-variation: 



Consider diatomic molecules with the unpaired elec- 
tron and the 2 £ ground state. Examples of such 
molecules include LaS, LaO, LuS, LuO, and YbF 
|73j . The hyperfine interval Ahf s is proportional to 
a 2 ZF rc i(aZ)fig nuc , where F rel is an additional relativistic 
(Casimir) factor |75j . The rotational interval A lot cx fi 
is roughly independent on a. If we find a molecule with 
Ahfs ~ A ro t, the splitting u> between hyperfine and rota- 
tional levels will depend on the following combination 



/x [a 2 F Te \(aZ) g nuc - const] 



Relative variation is then given by 



8lo 



Mils 



(2 + K) 



5a 5g n 



9v 



(39) 



(40) 



where the factor K comes from variation of F ve \(aZ), 
and for Z ~ 50, K w 1. As long as Ahf s /a; 1, we can 
neglect the last term in Eq. ( 40 1 . 



The data on hyperfine structure of diatomics are sparse 
and usually not very accurate. This hampers the search 
for molecules with strong cancelation of the types, dis- 
cussed here. Using data from [73] one can find that 
u) = (0.002 ± 0.01) cm" 1 for 139 La 32 S 05]. Note that 
for uj = 0.002 cm -1 the relative frequency shift is: 



Stu 
u> 



600 



5a 



(41) 



fi/(j,= (3.4 ±6.5) x lO^yr -1 . 



(38) 



This limit is significantly weaker than astrophysical limit 



With new data on molecular hyperfine constants appear- 
ing regularly, it is likely that other molecular candidates 
for such experiments will appear soon. 
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B. Molecules with cancelation between 
fine-structure and vibrational intervals 

The fine-structure interval u>f rapidly grows with the 
nuclear charge Z: 



uif ~ Z a , 



(42) 



On the contrary, the vibration energy quantum decreases 
with the atomic mass: 



^vib ~ M- 1 ' 2 ^ 1 ' 2 



(43) 



where the reduced mass for the molecular vibration is 
M r m p . Therefore, we obtain an equation Z = Z(M r ,v) 
for the lines on the plane Z,M r , where we can expect 
approximate cancelation between the fine-structure and 
vibrational intervals: 



ui = ujj — wcjvib w , v—1,2, 



(44) 



Using Eqs. ( 42 44 1 it is easy to find the dependence of 
the transition frequency on the fundamental constants: 



Slu 1 / 5a v 5\i 

— = — 2w f 1- -uj vih — 

oj oj \ a 2 



K 2 



5a 1 5a 
\ 

2 jJL 



(45) 



where the enhancement factor K — determines the 

UJ 

relative frequency shift for the given change of funda- 
mental constants. Large values of the factor K hint at 
potentially favorable cases for performing an experiment 
because it is usually preferable to have larger relative 
shifts. However, there is no strict rule that larger K is 
always better. In some cases, such as very close levels, 
this factor may become irrelevant [42]. Thus, it is also 
important to consider the absolute values of the shifts 
and compare them to the linewidths of the correspond- 
ing transitions. 

Because the number of molecules is finite we can not 
have lu — exactly. However, for many molecules we 
do have ui/uif <C 1 and \K\ ^> 1. Moreover, an addi- 
tional "fine tuning" may be achieved by selection of iso- 
topes and rotational, fi-doublet, and hyperfine compo- 
nents. Therefore, we have two large manifolds, the first 
one is built on the electron fine-structure excited state, 
and the second one is built on the vibrational excited 
state. If these manifolds overlap, one may select two or 
more transitions with different signs of lu. In this case ex- 
pected sign of the |w|-variation must be different (since 
the variation 5u> has the same sign) and one can elimi- 
nate some systematic effects. Such control of systematic 
effects was used in Refs. H3] [73] for transitions be- 
tween close levels in two dysprosium isotopes. The sign 
of energy difference between two levels belonging to dif- 
ferent electron configurations was opposite for the 163 Dy 
and 162 Dy transitions used in that work. 

In Table[l]we present the list of molecules from Ref. [73] 
where the ground state is split in two fine-structure levels 



TABLE I: Diatomic molecules with quasidegeneracy between 
the ground-state vibrational and fine-structure excitations. 
All frequencies are in cm -1 . The data are taken from Ref. [73] . 



Enhancement factor K is estimated using Eq. ( 45 1 



Molecule Electronic states 




W v ib 


K 


CI2 2 Il3/2.1/2 


645 


645.6 


1600 


CuS 2 n 


433.4 


415 


24 


IrC 2 A 5/2! 3/2 


3200 


1060 


160 


SiBr 2 n 1/2i3/2 


423.1 


424.3 


350 



and Eq. (44 1 is approximately fulfilled. The molecules 



ClJ and SiBr are particularly interesting. For both of 



them the frequency uj defined by Eq. (44 1 is of the or- 
der of 1 cm' 1 and comparable to the rotational constant 
B. This means that uj can be reduced further by proper 
choice of isotopes, rotational quantum number J and hy- 
perfine components. New dedicated measurements are 
needed to determine exact values of the transition fre- 
quencies and to find the best transitions. However, it 
is easy to find the necessary accuracy of the frequency- 
shift measurements. According to Eq. (45), the expected 
frequency shift is 



5uj = 2lo 



5 a 
a 



1 5/i 

4 jJL 



(46) 



Assuming 5a /a ~ 10' 15 and ujf ~ 500 cm" 1 , we obtain 
5ui ~ 10" 12 cm -1 ~ 3x 10 -2 Hz. In order to obtain simi- 
lar sensitivity comparing hyperfine transition frequencies 
for Cs and Rb one has to measure the shift ~ 10 -5 Hz. 



C. Molecular ion HfF+ 

The list of molecules in Table [I] is incomplete because 
of the lack of data in Ref. [73]. Let us briefly discuss 
one interesting case, which has been brought to atten- 
tion quite recently. The HfF + ion and other similar ions 
are being considered by E. Cornell's group at JILA for 
an experiment to search for the electric dipole moment 
(EDM) of the electron [7S] [77]. In this experiment, the 
ions are to be trapped in a quadrupole RF trap to achieve 
long coherence times. A similar experimental setup can 
be used to study possible time-variation of fundamental 
constants. A recent calculation by Petrov et al. [73] sug- 
gests that the ground state of this ion is 1 S + and the first 
excited state 3 Ai lies only 1633 cm" 1 higher. The calcu- 
lated vibrational frequencies for these two states are 790 
and 746 cm" 1 , respectively. For these parameters the vi- 
brational level v — 3 of the ground state is only 10 cm" 1 
away from the v = 1 level of the 3 Ai state. Thus, instead 
of Eq. (44 1 we now have: 



3, ,(1) 



7, ,(0) 
2 W vib 



0, 



(47) 
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where superscripts and 1 correspond to the ground 
and excited electronic states. The electronic transition 
with frequency lu c \ is not a fine-structure transition, and 
Eq. (42) is not applicable. Instead, by analogy with 
Eq. (|6| we can write: 



Wel = Wel.O + 1 X ■ 



a /a 



(48) 



In order to calculate the g-factor for HfF + ion one 
needs to perform relativistic molecular calculations for 
several values of a, which has not been done yet. How- 
ever, it is possible to make an order of magnitude es- 
timate using atomic calculation for the Yb + ion [2"4"j . 
According to Rcf. [78], the 1 T,^ - 3 Ai transition, to a 
first approximation, corresponds to the 6s - 5d transi- 
tion in the hafnium ion. It is well known that valence 
s- and d-orbitals of heavy atoms have very different de- 
pendence on a: while the binding energy of s-electrons 
grows with a, the binding energy of rf-clcctrons decreases 
US! HQ 121] ■ For the same transition in the Yb+ ion 
the Ref. [22] gives q s d = 10000 cm -1 . Using this value as 
an estimate, we can write by analogy with Eq. (451: 



5u> « 20000 cm- 1 (5a/a + 0.045 fi/^i) 



2000— +80^1, (49) 



(50) 



Assuming 8a/ a ~ 10 15 we obtain 5uj ~ 0.6 Hz. 



D. Estimate of the natural widths of the 
quasidegenerate states 

As we mentioned above, it is important to compare 
frequency shifts caused by time-variation of constants to 
the linewidths of corresponding transitions. First let us 
estimate the natural width r„ of the vibrational level v. 



T v = 



4uA 



vib I 



she 3 



[v\D\v-D 



(51) 



To estimate the dipole matrix element we can write: 



D 



dD(R) 



OR 



(R - Ro) 



R=Ro 



Do 
Ro 



(R - R ) , (52) 



where Do is the dipole moment of the molecule for equi- 
librium internuclear distance i?o- Using the standard 
expression for the harmonic oscillator, (u|a;|w — 1) = 
(hv/2muj) 1 ^ 2 , we get: 



3c 3 M r m p i?Q 



For the homonuclear molecule Cl 2 Dq 



(53) 



and expression 



Fi 



(53 I turns to zero. For the SiBr molecule, Eq. (53 1 gives 
10~ 2 Hz, where we assume D 2 /R$ ~ 0.1 e . 
Now let us estimate the width of the upper state 



of the fine-structure doublet 2 n!/2,3/2- By analogy with 
Eq. (51) we can write: 



r / = ^Jl( 2n 3/ 2 |^il 2 n 1 /2)| 



(54) 



The dipole matrix element in this expression is written 
in the molecular frame and we have summed over final 
rotational states. This matrix element corresponds to 
a spin-flip and turns to zero in the non-relativistic ap- 
proximation. Spin-orbit interaction mixes the 2 n 1 / 2 and 
2 Xi/ 2 states: 



2 n 1/2 ) _> | 2 n 1/2 ) + e| 2 Si/ 2 ) , 



and the matrix element in Eq. ( 54 ) becomes 79J : 



^ns/alDifni/^MCOlpilE) 



a 2 Z 2 
10(£ n - E s ) 



(55) 



(56) 



where Es is the energy of the lowest S-state. Substi- 
tuting Eq. ( 56 1 into Eq. ( 54 1 and using energies from 



Ref. [73], we get the following estimate for the molecules 

(57) 



Clt and SiBr 



F/ 



10" 2 Hz . 



Here we took into account that the unpaired electron 
in SiBr molecule is predominantly located at Si (Z=14) 
rather then at Br (Z=35). Because of this, the fine- 
structure splitting in SiBr is smaller than that of Cl^, 
where Z = 17 (see Table [I). 

We conclude that natural widths of the molecular lev- 
els considered here are of the order of 10 -2 Hz. This can 
be compared, for example, to the natural width 12 Hz 
of the level 2 D 5 / 2 of Hg + ion, which was used in atomic 
experiment [4 . For such narrow levels the lifetime may 
depend on the interaction with the black body radiation 
[5U] , According to this reference, the lifetimes of the ro- 
vibrational levels of polar molecules at room temperature 
vary from 1 s to 100 s. 



IX. EXPERIMENTS WITH Cs 2 AND Sr 2 

In this section we discuss two recently proposed ex- 
periments with cold diatomic molecules. The first one 
with Cs 2 molecule was proposed at Yale 44. 81 and the 
second experiment with Sr 2 molecule is in preparation at 
JILA [52]. 

The Yale experiment is based on the idea [35 to match 
an electronic energy with a large number of vibrational 



quanta. The difference with Eqs. (42 - 44 1 is that here 



electronic transition is between the ground state ^^j" and 
the 3 S+ state and, to a first approximation, its frequency 
is independent of a. The energy of this transition is about 
3300 cm -1 and the number of vibrational quanta needed 
to match this interval is on the order 100 (see Fig. [TJ . For 
the vibrational quantum number v ~ 100 the density of 
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Internuclear Separation in a.u. 



FIG. 1: Levels 3 E+ and X E+ in Cs 2 molecule (figure from Ref. [83]). 



levels is high due to unharmonicity and it is possible to 
find very close levels of two different potential curves. 
This leads to enhanced sensitivity to variation of p, as in 
Eq. (44 1. Cold CS2 molecules can be produced in a par- 



ticular quantum state by photoassociation of Cs atoms 
in a trap. 



Let us estimate the sensitivity of this experiment to 
variation of a and p. For the electronic transition energy 
we can use Eq. ( 48 1 . If we neglect unharmonicity, we can 



write the transition frequency between close vibrational 
levels of the two electronic terms in the form 

0J = ojd.o + qx+ (v 2 + \) w V ib,2 - {v\ + \) Wvib,i) (58) 

where v-i <C v\. The dependence of this frequency on 
constants is given by: 



8a Wei o 8p 

2q ^ 

a 2 p 



(59) 



where we took into account that u> <C w e i,o- A very rough 
estimate of the factor q can be done in the following way. 
For the ground state of atomic Cs the g-factor is about 
1100 cm -1 , which is close to \a 2 Z 2 e§ s , where eg s is the 
ground-state binding energy. If we assume that the same 
relation holds for the electronic transition in molecule, we 



get |g| 



Wei,o ~ 120 cm . Using this estimate 



and Eq. (59 1 we get: 
8u> ~ 



-240 



5 a 



1600 



8p 
A 4 



(60) 



where we assume that relativistic corrections reduce dis- 
sociation energy of the molecule, so q is negative. This 
estimate shows that the experiment with Cs2 is mostly 
sensitive to variation of /x. 



Estimate (60 1 is obtained in the harmonic approxima- 
tion. As mentioned above, for high vibrational states 
real potential is highly unharmonic. This significantly 
decreases the sensitivity of this experiment compared to 



the naive estimate ( 60 1 . It can be easily seen cither from 



the WKB approximation [3U [5T] , or from an analytical 
solution for the Morse potential [52] • Quantization condi- 
tion for vibrational spectrum in the WKB approximation 
reads: 



rt 2 



Ri 



y/2M(U(r)-E n )dr = (v + §) tt . 



Differentiating this expression in [i we get: 



6E V 



2p(E v ) (i ' 



(61) 



(62) 



where p(E v ) = (dE v /dv) 1 w (E v — E v _ x ) 1 is the level 
density. For the harmonic part of the potential, p = const 
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and the shift 6E V grows linearly with v, but for vibra- 
tional states near the dissociation limit the level density 
p(E) — ► oo and SE V — ► 0. Consequently, maximum 
sensitivity ~ 1000 cm -1 is reached at v s» 60, and rapidly 
drops for higher v. At present the group at Yale has 
found a conveniently close vibrational level of the upper 
3 E„ state for v = 138, but where the sensitivity is only 
~ 200 cm -1 [ST]. There are still good chances that there 
are other close levels with smaller v, where the sensitivity 
may be several times higher. 

It is important that because of unharmonicity, the sen- 
sitivity to variation of a also decreases compared to the 
estimate (60 1. The reason for this is the following. For 



the highest vibrational levels of the ground state, as well 
as for all levels of the upper (weakly bound) state, the 
separation between nuclei is large, R > 12 a.u. (sec 
Fig. [IJ. Thus, both electronic wave functions are close 
to either symmetric (for ^jj") or antisymmetric combi- 
nation (for 3 S^) of atomic 6s functions: 

^ g , u (r 1} r 2 ) » i (6s a (r 1 )6s i '(r 2 ) ± 6s b ( ri )6s a (r 2 )) . 

(63) 

Therefore, all rclativistic corrections are (almost) the 
same for both states. 

Similar conclusions can be reached from the analysis 
of the Morse potential: 



U AI (r)=d(l-e- a( - r - ro ^y -d. 



(64) 



The eigenvalues for this potential are given by the ana- 
lytical expression: 



E v = uj (v + |) - 



-d, 



(65) 



where ujq — 2nay/2d/M and the last eigenvalue En is 
found from the conditions En+i < En and En-i < En- 
Obviously, En is very close to zero and is practically in- 
dependent from any parameters of the model. Therefore, 
it is also insensitive to variation of constants. 

We see that highest absolute sensitivity is reached for 
vibrational levels somewhere in the middle of the poten- 
tial curve. However, in this part of the spectrum there are 
no close levels of different nature to maximize the relative 
sensitivity 5w/uj. One can still use frequency combs to 
perform high- accuracy measurements. This idea is used 
in the resent proposal by Zelevinsky et al. (author?) 
[52"] , who suggest to use an optical lattice to trap Sr2. 
These molecules are formed by photoassociation in one 
of the uppermost vibrational levels of the ground elec- 
tronic states (see Fig. [2]). As we saw above, this level is 
not sensitive to the variation of /i. At the next stage, a 
Raman transition is driven to one of the most sensitive 
levels in the middle of the potential well. This way it is 
possible to get the highest possible absolute sensitivity 
for a given molecule. Unfortunately, the dissociation en- 



ergy for Sr2 is only about 1000 cm -1 , which is 3 times 
smaller than for Cs2. Because of this, the highest sensi- 
tivity for the Sr 2 molecule is about 270 cm -1 , i.e. only 
slightly higher than the sensitivity of the v = 138 level 
in Cs2- Therefore, it may be useful to try to apply this 
scheme to some other molecule with larger dissociation 
energy. Note that in the experiment with Sr 2 , the sensi- 
tivity to a- variation is additionally suppressed by a factor 
(38/55) 2 w 1/2 because of the smaller Z. 



X. EXPERIMENTS WITH HYDROGEN 
MOLECULAR IONS H+ AND HD+ 

Hydrogen molecular ions are very attractive for funda- 
mental studies because of their theoretical simplicity and 
experimental possibility of their cooling and trapping. 
Using H+ and HD+ ions for studying time-variation 
of electron-to-proton and proton-to-deuteron mass ra- 
tios fi — m e /m p and m p /md has been suggested in 
Refs. [HH [55] ■ Because of the unharmonicity, the ratio 
of the two vibrational transitions with very different vi- 
brational quantum numbers is /x-dependent [53]. There 
is no enhancement of the relative effect here, but the 
lines are very narrow and high-precision measurements 
are possible using frequency combs. 

Recently HD + ion has been cooled to 50 mK and 
trapped in a linear rf trap [86J . This allowed to measure 
the rovibrational transition v,N = 0, 2 — > v',N' = 4,3 
with an absolute accuracy of 0.5 MHz. Using sensitiv- 
ity coefficient from [85] one can see, that this accuracy 
translates in to 5 x 10~ 9 (5 ppb) accuracy for /i. Note 
that modern molecular theory of HD + has comparable 
accuracy |87j . Thus, a direct comparison between theory 
and experiment allows to determine the absolute value of 
/i to 5 ppb. 



XL CONCLUSIONS 

We have seen that both diatomic and polyatomic 
molecules are used in astrophysics to study possible vari- 
ation of the electron-to-proton mass ratio \i on a time 
scale from 6 to 12 billion years. Results of these stud- 
ies are inconclusive, see Eqs. (15 1, (19 1, and (36). The 



situation is similar for the astrophysical search for co- 
variation. In principle, all these results can be explained 
by complex evolution of fi and a in space and time. Or, 
more likely, there are some systematic errors, which are 
not fully understood. Therefore it is extremely impor- 
tant to supplement astrophysical studies with laboratory 
measurements of present-day variation of these constants. 
This work is currently going on in many groups. Most of 
them use atomic frequency standards and atomic clocks. 
In this chapter we discussed several resent ideas and pro- 
posals on how to increase the sensitivity of laboratory 
tests by using molecules instead of atoms. 

The only molecular experiment |47l 148] . which has 
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FIG. 2: The scheme for Raman spectroscopy of Sr2 ground-state vibrational spacings. A two-color photoassociation pulse 
prepares molecules in the v = Umax — 2 vibrational state (denoted on the plot as v = —3). Subsequently, a Raman pulse couples 
the v — —3 and v — 27 states via v' « 40 level of the excited 0^ state (figure from Ref. [82]). 



reached the stage of placing the limit on the time- 
variation of fundamental constants (37 1, used a super- 



sonic molecular beam of SF 6 . Even though this experi- 
ment is less sensitive than the best atomic experiments, 
it constrains a different combination of fundamental con- 
stants. This allows to combine it with the results of 
atomic-clock experiments [U [35J [37] to place the most 
stringent laboratory limit (38 1 on time- variation of /i. 
The linewidth in this experiment, T w 200 Hz, was de- 
termined by the time-of-flight through the 1-m Ramsey 
interferometer. A similar problem with the linewidth has 
prevented the use of the ND3 beam to perform com- 
petitive experiment on time- variation [46] , Using cold 
molecules would allow to reduce the linewidth by several 
orders of magnitude and drastically raise the sensitivity 
of molecular experiments. 

We have seen that for such diatomic radicals as Clj" 
and SiBr there are narrow levels of different nature sep- 
arated by intervals < 1 cm . The natural widths of 



these levels are on the order of 10~ 2 Hz. This is compa- 
rable to the accuracy necessary to reach the sensitivity 
of 5a/a ~ 10~ 15 , similar to that of the best modern lab- 
oratory tests. In the high-precision frequency measure- 
ments, the measurement accuracy is typically few orders 
of magnitude better than the linewidth. Of course, in 
order to benefit from such narrow lines, it is crucial to be 
able to cool and trap the molecules. In this respect the 
ion CI 2 looks more promising. 

Even higher sensitivity to the temporal variation of a 
can be found in HfF + and similar molecular ions, which 
are being considered for the search of the electron EDM 
at JILA [76l [771 EH] ■ The transition amplitude between 
and 1 S of HfF + ion is also suppressed. The transi- 
tion width is larger than for ClJ and SiBr because of the 
larger value of Z and higher frequency uif. In Ref. [78] 
the width of 3 Ai state was estimated to be about 2 Hz. 
This width is also of the same order of magnitude as the 
expected frequency shift for 5a /a ~ 10~ 15 . At present 
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not much is known about these molecular ions. More 
spectroscopic and theoretical data are needed to estimate 
the sensitivity to a-variation reliably. We hope that this 
review may stimulate further studies in this direction. 
Additional advantage here is the possibility to measure 
electron EDM and a-variation using the same molecule 
and a similar experimental setup. 

Preliminary spectroscopic experiment with the CS2 
molecule has been recently finished at Yale |81| . The 
electron transition in CS2 goes between the 3 E+ and 
states, and, to a first approximation, is independent of 
a. On the other hand the sensitivity to /i may be en- 
hanced because of the large number of vibrational quanta 
needed to match the electronic transition. However, the 
unharmonicity of the potential curve near the dissocia- 
tion limit suppresses this enhancement for very high vi- 
brational levels. As a result, the sensitivity to variation 
of fx for the v — 138 level is about the same as in Eq. (46 1. 



It is possible that there are other close levels with smaller 
vibrational quantum number v and, consequently, with 
higher sensitivity. Even if such levels are not found, the 
experiment with the v = 138 level may improve present 
limit on variation of /i by several orders of magnitude. 

An experiment with the Sr2 molecule was recently pro- 
posed at JILA [82] . This experiment potentially has sim- 



ilar sensitivity to variation of /i as the experiment with 
Cs2 and both of them are complementary to the experi- 
ments with molecular radicals, which are mostly sensitive 
to a-variation [49 . 

Finally, we have seen that the inversion spectra of such 
polyatomic molecules as NH 3 and ND 3 are potentially 
even more sensitive to variation of fi. This has already 
been used in astrophysics to place the most stringent 
limit ( 36 ) on the time- variation of /j, on the cosmologi- 



cal timescale. Corresponding laboratory experiments re- 
quire very slow molecular beams, fountains, or molecular 
traps. The work in this direction is going on |46j . 

To conclude this chapter, we see that this field is 
rapidly developing and new interesting results can be ex- 
pected in the near future. 
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